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Abstract 

We introduce a microscopic model which describes the dynamics of each dealer 
in multiple foreign exchange markets, taking account of the triangular arbitrage 
transaction. The model reproduces the interaction among the markets well. We 
explore the relation between the parameters of the present microscopic model and 
the spring constant of a macroscopic model that we proposed previously. 
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1 Introduction 



Analyzing correlation in financial time series is a topic of considerable in- 
terest [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. We recently pointed out 
[19,20,21] that, in the foreign exchange market, a correlation among exchange 
rates can be generated by triangular arbitrage transactions. The triangular 
arbitrage transaction is a financial activity that takes advantage of the three 
foreign exchange rates among three currencies [19,22,23]. It makes the product 
of the three foreign exchange rates converge to its average, thereby generating 
an interaction among the rates. 

In order to study effects of the triangular arbitrage on the fluctuations of 
the exchange rates, we introduced [19] a stochastic model describing the time 
evolution of the exchange rates with an interaction. The model successfully 
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described the fluctuation of the data of the real market. The model is phe- 
nomenological; i.e. it treats the fluctuations of the rates as fluctuating particles 
and the interaction among the rates as a spring. We refer to this model as the 
'macroscopic model' hereafter. 

The purpose of this paper is to understand microscopically effects of the tri- 
angular arbitrage on the foreign exchange market. For the purpose, we intro- 
duce a new model which focuses on each dealer in the markets; we refer to 
the new model as the 'microscopic model' hereafter. We then show the rela- 
tion between the macroscopic model and the microscopic model through an 
interaction strength which is regarded as a spring constant. 

The paper is organized as follows. We first review the macroscopic model 
of the triangular arbitrage in Section 2. Second, in Section 3, we introduce 
the microscopic model which focuses on the dynamics of each dealer in the 
markets. The model reproduces the interactions among the markets well. We 
explore the relation between the spring constant of the macroscopic model 
and the parameters in the microscopic model in Section 4. We summarize the 
paper in Section 5. 



2 Review of macroscopic model of triangular arbitrage 



The triangular arbitrage is a financial activity that takes advantage of three 

exchange rates. When a trader exchanges one Japanese yen to some amount 
of US dollar, exchanges the amount of US dollar to some amount of euro and 
exchanges the amount of euro back to Japanese yen instantly at time t, the 
final amount of Japanese yen is given by 

l^=llr,it), (1) 

x=l 

where 

"■'^^^^ yen-dollar ask (t) 

r2it) = ^-n r-TT (3) 

dollar-euro ask [t) 

r3(t) = yen-euro bid (t). (4) 

Here, 'bid' and 'ask,' respectively, represent the best bidding prices to buy 
and to sell in each market. If the rate product fi is greater than unity, the 
trader can make profit through the above transaction. This is the triangular 
arbitrage transaction. Once there is a triangular arbitrage opportunity, many 
traders will make the transaction. This makes /i converge to a value less than 
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(a) (b) 

Fig. 1. The probability density function of the rate product /i [19]. (b) is a 
semi- logarithmic plot of (a). The shaded area represents triangular arbitrage op- 
portunities. The data were taken from January 25 1999 to March 12 1999. 

unity, thereby eliminating the opportunity. Triangular arbitrage opportunities 
nevertheless appear, because each rate n fluctuates strongly. 

The probability density function of the rate product has a sharp peak and 
fat tails (Fig. 1). It means that the fluctuations of the exchange rates have 
correlation that makes the rate product converge to its average (//) ~ 0.99998. 
The average is less than unity because of the spread; the spread is the difference 
between the ask and the bid prices and is usually of the order of 0.05% of the 
price. 

For later convenience, we here define the logarithm rate product u as the 
logarithm of the product of the three rates: 

3 3 

u{t) ^lnl[r,{t) = Y.\nr,{t). (5) 

x=l x=l 

There is a triangular arbitrage opportunity whenever this value is positive. 

We can define another logarithm rate product u', which has the opposite 
direction of the arbitrage transaction to u, that is, from Japanese yen to euro 
to US dollar back to Japanese yen: 



u\t) = j2\nr\{t), (6) 

x=l 

where 

r'i{t) = yen-doUar bid (t) (7) 
r'2{t) = doUar-euro bid (t) (8) 

- yen-eJaskW 

This logarithm rate product v' will appear in Section 3.2. 
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Fig. 2. The probability density 
function of u [19]. The circle (o) 
denotes the real data and the 
solid line denotes our simulation 
data. We fix the time step at 
T = 60 [sec] and hence use the 
spring constant k = 0.17 for our 
simulation. The simulation data 
fit the real data well. 

In one of our previous works [19], we constructed a stochastic model of the 
time evolution of foreign exchange rates that takes account of the effect of the 
triangular arbitrage transaction. The basic equation of the model is the time 
evolution of the logarithm of each rate: 



A schematic image of the 
[24]. The three random 
walkers with the restoring force 
working on the center of gravity. 



lnr^(t + T) = Inr^(t) + ri^{t) + g{v{t)), 



[X 



1,2,3) 



(10) 



where T is a time step which controls the time scale of the model; we later use 
the actual financial data every T[sec]. The variable 7],^. denotes an independent 
fluctuation that obeys a truncated Levy distribution [9,6,25] and g represents 
an interaction function defined by 



fill 



where /c is a positive constant which specifies the interaction strength and {v) 
is the time average of v. The time-evolution equation of the logarithm rate 
product v is given by summing Eq. (10) over all x: 



v{t + T)- {u) = (1 - 3k){u{t) - ^)) + E vAt)- 



[12) 



x=l 



The model equation (12) well describes a fat-tail probability distribution of 
u{t) of the actual market (Fig. 2) [19] as well as a negative auto-correlation of 
the price fiuctuation [20]. 

From a physical viewpoint, we can regard [24] the model equation (10) as 
a one-dimensional random walk of three particles with a restoring force, by 
interpreting Inr^; as the position of each particle (Fig. 3). The logarithm rate 
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Fig. 4. The spring constant k as a function of the time step T. The spring constant 
k increases with the time step T. 

product is the summation of InVx, hence is proportional to the center of 
gravity of the three particles. The restoring force glu) makes the center of 
gravity converge to a certain point (u). The form of the restoring force (11) is 
the same as that of the harmonic oscillator. Hence we can regard the coefficient 
k as a spring constant. 

The spring constant k is related to the auto-correlation function of u as follows 
[19]: 

i-3k- cm = Mt + T).(t))-Mt))^ .^3. 

Using Eq. (13), we can estimate the spring constant k{T) from the real data 
series as a function of the time step T (Fig. 4). The spring constant k increases 
with the time step T. We fixed the time step at T = 60 [sec] and hence used 
the spring constant k — 0.17 for our simulation. We will come back to this 
point later in Section 4. 



3 Microscopic model of triangular arbitrage 



We here introduce a microscopic model which describes interactions among 
foreign exchange markets. The model focuses on the dynamics of each dealer 
in the market. 

In order to describe each foreign exchange market microscopically, we use Sato 
and Takayasu's dealer model (the ST model) [26] , which reproduces the power- 
law behavior of price changes in a single market well. Although we focus on the 
interactions among three currencies, two of the three markets can be regarded 
as one effective market [21]; i.e. the yen-euro rate and the euro-dollar rate are 
combined to an effective yen-dollar rate. In terms of the macroscopic model, we 
can redefine a variable r2 as the product of r2 and r^. Then the renormalized 
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Fig. 5. A schematic image of a transaction of the ST model. Only the best bidders are 
illustrated in order to simplify the image. The circles denote the dealers' bidding 
price to buy and the squares denote the dealers' bidding price to sell. The filled 
circles denote the best bidding price to buy, max{i?j}, and the gray squares denote 
the best bidding price to sell, min{Sj} + A. In (a), the condition (15) is not satisfied, 
and the dealers, following Eq. (17), change their relative positions by Oj. Note that 
the term cAP does not depend on i; hence it does not change the relative positions 
of dealers but change the whole dealers' positions. In (b), the best bidders satisfy 
the condition (15). The price P is renewed according to Eq. (16), and the buyer 
and the seller, respectively, become a seller and a buyer according to Eq. (18). 

variable r2 follows a similar time-evolution equation. We therefore describe 
triangular arbitrage opportunities with only two interacting ST models. 



3.1 Sato and Takayasu's dealer model (ST model) 



We first review the ST model briefly [26] (Fig. 5). The basic assumption of 
the ST model is that dealers want to buy stocks or currencies at a lower price 
and to sell them at a higher price. There are N dealers; the ith dealer has 
bidding prices to buy, Bi{t), and to sell, Si{t), at time t. Let us assume that 
the difference between the buying price and the selling price is a constant 
A = Si{t) — Bi{t) > for all i, in order to simplify the model. 

The model assumes that a trade takes place between the dealer who proposes 
the maximum buying price and the one who proposes the minimum selling 
price. A transaction thus takes place when the condition 

max{5i(t)} > mm{Si{t)} (14) 

or 

max{Bi{t)} - mm{Bi{t)} > A (15) 

is satisfied, where max{-} and min{-}, respectively, denote the maximum and 
the minimum values in the set of the dealers' buying threshold {Bi{t)}. The 
logarithm of the rates in the actual market, —Inr^ and Inr^., correspond to 
max{i?i} and minlS*,}, respectively. 
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The market price P{t) is defined by the mean vahie of max{Bi} and minlS'j} 
when the trade takes place. The price P{t) maintains its previous value when 
the condition (15) is not satisfied: 



P{t) 



(max{5j(i)} + min{Si{t)})/2, if the condition (15) is satisfied, 
P{t — 1), otherwise. 



(16) 



The dealers change their prices in a unit time by the following deterministic 



where aj(t) denotes the ith dealer's characteristic movement in the price at 
time t, AP(t) is the difference between the price at time t and the price at 
the time when the previous trade was done, and c(> 0) is a constant which 
specifies dealers' response to the market price change, and is common to all 
of the dealers in the market. The absolute value of a dealer's characteristic 
movement ai(t) is given by a uniform random number in the range [0, a) and 
is fixed throughout the time. The sign of aj is positive when the ith dealer is 
a buyer and is negative when the dealer is a seller. The buyer (seller) dealers 
move their prices up (down) until the condition (15) is satisfied. Once the 
transaction takes place, the buyer of the transaction becomes a seller and the 
seller of the transaction becomes a buyer; in other words, the buyer dealer 
changes the sign of from positive to negative and the seller dealer changes 
it from negative to positive: 



The initial values of {Bi} are given by uniform random numbers in the range 
(—A, A). We thus simulate this model specifying the following four parameters: 
the number of dealers, A^; the spread between the buying price and the selling 
price. A; the dealers' response to the market price change, c; and the average 
of dealers' characteristic movements in a unit time, a. 

The ST model well reproduces the power-law behavior of the price change 
when the dealers' response to the market change c > c*, where c* is a critical 
value to the power-law behavior. The critical point depends on the other 
parameters; e.g. c* ~ 0.25 for N = 100, A = 1.0 and a = 0.01 [26]. 



rule: 



Bi{t + 1) = Bi{t) + Uiit) + cAP{t), 



(17) 




—ai{t) for the buyer and the seller, 
ai{t) for other dealers. 



(18) 
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3.2 Microscopic model of triangular arbitrage: interacting two systems of the 
ST model 



We now describe our microscopic model as a set of the ST models. In order 

to reproduce effects of the triangular arbitrage, we prepare two systems of the 
ST model, the market X and the market Y. As is noted above, we prepare 
only two markets to reproduce the effect of the triangular arbitrage because 
we regard two of the three markets as one effective market. 

The dealers in the markets X and Y change their bidding prices according to 
the ST model as follows: 



where X and Y denote the markets X and Y, respectively. An intra-market 
transaction takes place when the condition 



is satisfied. We assume that A is common to the two markets. The price Px{t) 
is renewed in analog to the ST model: 



where x — X or Y. 

We here add a new mter-market transaction rule which makes the systems 
interact. The arbitrage transaction can take place when one of the conditions 



is satisfied (see Fig. 6). When the conditions (21) and (23) or (24) are both 
satisfied simultaneously, the condition (21) precedes. 

Note that the arbitrage conditions z^x > and z/y > in the microscopic model 
correspond to the arbitrage condition z/ > in the actual market, where u is 
defined by Eq. (5). We assume that the dealers' bidding prices {Bi} and {Si} 
correspond to the logarithm of the exchange rate, Inrj. Therefore, maxjSj x} 
may be equivalent to — In (yen-dollar ask) while minj^j^y} may be equivalent 
to ln(dollar-euro ask) — ln(yen-euro bid), and hence ux niay be equivalent to 
I/. More precisely, the direction of the arbitrage transaction determines which 
of the quantities, ux or uy, corresponds to the logarithm rate product u. There 



Bi,x{t + 1) = Bi^xit) + ai,x{t) + cAPxit) and 
Bixit + 1) = Bi,Y(t) + ai,Y(t) + cAPrit), 



(19) 
(20) 



max{i?i,^(t)} > mm{Si^x{t)} , x ^ X or Y 



(21) 




ux = max{Bi,x{t)} - (min{5i,y(t)} + A) > 
uy = max{A,y(i)} - (min{Si,x(i)} + A) > 



(23) 
(24) 
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Fig. 6. A schematic image of the transactions. Only the best bidders in the markets 
are illustrated in order to simplify the image. The circles and the squares denote 
the dealers' bidding price to buy and to sell. The filled circles denote the best 
bidding prices to buy in the markets, max{5j_x} and max{f?j y}, and the gray 
squares denote the best bidding prices to sell in the markets, mm{Bi^x} + A and 
mm{Bi Y} + A. In the case (a), any of the conditions (21), (23) and (24) are not 
satisfied. The buyers move their prices up, and the sellers move their prices down. 
In the case (b), the dealers in the market X satisfy the condition (21); hence the 
intra-market transaction takes place. The price in the market X, Px, is renewed, 
and the buyer and the seller of the transaction become a seller and a buyer, respec- 
tively. In the case (c), the seller in the market X and the buyer in the market Y 
satisfy the condition (24); hence the arbitrage transaction takes place. The price 
Px in the market X becomes m.m{Bi^x} + A, and the price Py in the market Y 
becomes maxjSj^y}. The buyer and the seller of the transaction become a seller and 
a bTiycr, respectively. The arbitrage transaction thus makes the interaction between 
the markets X and Y. 

are two directions of the triangular arbitrage transaction. The definition (5) 
specifically has the direction of Japanese yen to US dollar to euro to Japanese 
yen. As is mentioned in Section 2, we can define another logarithm rate product 
u' in the actual market which has the opposite direction to u, Japanese yen 
to euro to US dollar to Japanese yen. Hence, if the logarithm rate product u 
in the actual market corresponds to vx in Eq. (23), u' corresponds to —i/y in 
Eq. (24). 

The procedures of the simulation of the microscopic model are as follows 
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(a) 



(b) 



Fig. 7. The distributions of vx and i/y- The parameters are fixed to N = 100, 
a = 0.01, A = 1.0 and (a) c = 0.0, (b) c = 0.3, and are common to the market 
X and the market Y. The sohd hne denotes vx and the dashed hne denotes uy in 
each graph. 



(1) Prepare two systems of the ST model, the market X and the market 
y, as described in Section 3.1. The parameters are common to the two 
systems. 

(2) Check the condition (21) and renew the prices by Eq. (22). If the condition 
(21) is satisfied, skip the step 3 and proceed to the step 4. Otherwise, 
proceed to the step 3. 

(3) Check the arbitrage conditions (23) and (24). If the condition (23) is satis- 
fied, renew the prices Px{t) and -Py(t) to max{5j^x(^)} and mm{Biy{t)}+ 
A, respectively. If the condition (24) is satisfied, renew the prices Px{t) 
and -Py(t) to min{Bi x{t)} + A and max{i?j y (t)}, respectively. If both of 
the conditions in (23) and (24) are satisfied, choose one of them with the 
probability of 50% and carry out the arbitrage transaction as described 
just above. If the arbitrage transaction takes place, proceed to the step 
4; otherwise skip the step 4 and proceed to the step 5. 

(4) Calculate the difference between the new prices and the previous prices, 
APx{t) = Pxit) - Pxit - 1) and APy(t) = Py(t) - Pyit - 1), and use 
them in the Eqs. (19) and (20), respectively. Change the buyer and the 
seller of the transaction to a seller and a buyer, respectively. In other 
words, change the signs of Oj^x and aiy of the dealers who transacted. 

(5) If any of the conditions (21), (23) and (24) are not satisfied, maintain 
the previous prices, Pxii) = Pxii — 1) and -Py(t) = Py(t — 1), as well 
as the previous price differences, APx(^) = ^Px{t — 1) and APyit) — 
APrit-l). 

(6) Change the dealers' bidding prices following Eqs. (19) and (20). 

(7) Repeat the steps from 2 to 6. 

The quantities ux and uy are shown in Fig. 7. (Following Ref. [26], we set 



(Fig. 6): 
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(a) 



(b) 



Fig. 8. The probability distribution of the difference between the logarithm rate 
product and its average for (a) the actual data, — f peak I and for (b) the data of 
the microscopic model, \i'x — ^a;,peaic| for c — 0.3, where J^peak and I'xjpeak 

denote the 

peak position of respective data. The solid lines represent the part of the positive 
difference and the dotted lines represent the part of the negative difference, in the 
both graphs. We can see that the probability distribution of the logarithm rate prod- 
uct v has a skewness around its average, and the microscopic model qualitatively 
reproduces it well. 

the parameters common to the two markets X and Y: N — 100, a — 0.01 
and A = 1.0.) In Fig. 7(b) for c = 0.3, the fat-tail behavior of the price 
difference i^x is consistent with the actual data as well as with the macroscopic 
model in Fig. 2. Furthermore, Px reproduces the skewness of the actual data, 
which cannot be reproduced by the macroscopic model (Fig. 8). Note that the 
skewness of uy is consistent with the behavior of u'. 



4 The microscopic parameters and the macroscopic spring con- 



In this Section, we discuss the relation between the macroscopic model and 
the microscopic model through the interaction strength, or the spring constant 



In the microscopic model, we define the spring constant /Cmicro, which corre- 
sponds to the spring constant k of the macroscopic model, as follows: 



stant 



k. 



micro — „ 




{ux{t + i)Mt))-{Mt))' 
{Mm - {MtW 



) 



(25) 



Figure 9 shows the estimate (25) as a function of several parameters. 
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(a) 



(b) 



Fig. 9. The spring COnstclIlt ^micro ^ 

function of parameters. The panel (a) 
shows the dependence on N. The other parameters are fixed to A = 1.0 and 
a = 0.0001,0.001 and 0.01 for the circles, the squares and the diamonds, respec- 
tively, and c = 0.3 for all the plots. The panel (b) shows the dependence on a. The 
other parameters are fixed to A = 1.0 and N = 100, 1000 and 10000 for the circles, 
the squares and the diamonds, respectively, and c = 0.3 for all the plots. 

Remember that, in the macroscopic model, the spring constant k depends on 
the time step T (see Fig. 4). The spring constant of the microscopic model 
fcmicro also depends on a time scale as follows. The time scale of the ST model 
may be given by the following combination of the parameters [26]: 



(n) 



3A 



(26) 



where n denotes the interval between two consecutive trades. Hence, the in- 
verse of Eq. (26), 

/ = l/(n) ^ ^, (27) 



3A 



is the frequency of the trades. 



Although there are four parameters N, a, A and c, we change only three 
parameters N, a, and c and set A = 1.0, because only the ratios N/A and 
a/ A are relevant in this system. The ratio N/A controls the density of the 
dealers and a/A controls the speed of the dealers' motion on average. Hence, 
we set A = 1.0 and change the other parameters hereafter. 

We plot the spring constant /Cmicro as a function of the trade frequency / = 
Na/3A in Fig. 10. The plots show that the spring constant A;inicro(-^) A) can 
be scaled by the trade frequency / well. 

In order to determine a reasonable range of the parameters, let us consider the 
situation in Fig. 11, where the arbitrage transaction is about to take place. 
At the moment, the positions of the second best bidders (hexagons) in the 
markets X and Y are, on average, A/N away from the prices transacted, Px 
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(a) (b) 
Fig. 10. The scaling plot of the spring const 8,nt A^micro 

as a function of the trade 

frequency / = Na/3A. The vertical axes are displayed in the logarithmic scale. The 
dealers' response to the price change, c, is fixed to 0.0 in (a) and 0.3 in (b). We 
fix a = 0.0001,0.001 and 0.01 and change N (open circles, squares, and diamonds, 
respectively) and A'^ = 100, 1000 and 10000 and change a (crosses, filled circles and 
triangles, respectively), while A is fixed to 1. Note that all points collapse onto a 
single curve. The spring constant fcmicro is scaled by /, and decays exponentially in 
both of the plots (a) and (b). 



□ f a/2 I 




Fig. 11. A schematic image of the second best bidders' motion. The circles and the 
squares denote the dealers' bidding price to buy and to sell. The filled circles and 
the gray squares represent the best bidding prices to buy and sell, respectively. The 
hexagons denote the second best bidding prices to buy (the filled one) and to sell 
(the gray one). 
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Fig. 12. The dependence of /o(c) on c, estimated by fitting the data in Fig. 10 as 
well as the same plots for different values of c. 

and Py. In the next step, the second best bidders in the markets X and Y will 
move by a/2 on average toward to the prices Px and Py, respectively. The 
next transaction will be carried out probably by the second best bidders. For 
a/2 > A/N, the prices of the transactions may move away from each other. 
The arbitrage transaction cannot bind the two prices of the markets X and 
Y enough and the two prices Px and Py fluctuate rather freely. It is not a 
realistic situation. Therefore, the condition 



/ 



Na ^2 
3A - 3 



(2J 



should be satisfied for the real market to be reproduced. On the other side, 
the simulation data have too large errors in the region / < 1/3 because the 
transaction rarely occurs. We hence use the data in the region 1/3 < / < 2/3 
hereafter. 



The spring constant /Cmicm decays exponentially in the range 1/3 < / < 2/3 
in both of the plots (a) and (b) of Fig. 10, having different slopes. Hence we 
assume that the spring constant decays as 

^micro ^ ^ f/foi )^ (29) 

where /o(c) denotes the characteristic frequency dependent on c. The estimates 
of the characteristic frequency /o(c) are shown in Fig. 12 as a function of c. 
The characteristic frequency /o(c) thus estimated decays linearly with c. The 
reason why /o(c) behaves so is an open problem. 

In Fig. 13, we plot the same data as in Fig. 4, but by making the horizontal axis 

the trade frequency /real- In order to compare it with Fig. 10 quantitatively, we 
used the time scale Tj-eai = 7[sec]; the interval between two consecutive trades 
in the actual foreign exchange market is, on average, known to be about 7[sec] 
[27]. The spring constant in the actual market k decays exponentially with the 
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Fig. 13. We plotted the same data as in Fig. 4, but the horizontal axis here is the 
trade frequency scaled by the realistic time scale of the trades, Tj-eai = 7 [sec]. 



trade frequency /real, which is consistent with that of the microscopic model 
shown in Fig. 10. The real characteristic frequency in Fig. 13, however, is quite 
different from that of the microscopic model plotted in Fig. 10. This is also 
an open problem. 

We have so far discussed the microscopic model consisting of two ST models. 
The model well describes the actual data qualitatively. We here note that we 
can reproduce the actual data by preparing straightforwardly a microscopic 
model consisting of three ST models (see Fig. 14). 




V 
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(a) (b) 

Fig. 14. The distributions of the logarithm rate products ux and z^y for the mi- 
croscopic model consisting of three ST models. The solid lines represent vx and 
the dashed lines represent i/y- The parameters are common to the three markets as 
N = 100, A = 1.0, a = 0.01 and (a) c = 0.0, and (b) c = 0.1. In (b) the model weh 
reproduces the fat-tail behavior of the actual data. 
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5 Summciry 



We first showed in Ref. [19] that triangular arbitrage opportunities exist in the 
foreign exchange market. The rate product /i fluctuates around its average. 
We then introduced in Ref. [19] the macroscopic model including the interac- 
tion caused by the triangular arbitrage transaction. Inspite of its simpleness, 
the maacroscopic model reproduced the actual behavior of the logarithm rate 
product 1/ well. We finally introduced here the microscopic model, which con- 
sists of two systems of the ST model. The microscopic model also reproduced 
the actual behavior of the logarithm rate product u well. The microscopic 
model can describe more details than the macroscopic model, in particular, 
the skewness of the distribution of the logarithm rate product u. We then 
explored the relation between the spring constant of the macroscopic model 
and the parameters in the microscopic model. The spring constant of the mi- 
croscopic model /Cmicro can be scaled by the trade frequency /, and it decays 
exponentially with /, which is consistent with the spring constant of the actual 
market k. 
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